Astronomy & Astrophysics manuscript no. Bonfini'Paolo 


© ESO 2009 


July 28, 2009 





ON 

o 
o 

(N 



00 
(N 

6 



Spectroscopy of the bright optical counterparts of X-ray sources 

in the direction of M 31. II 
(Research Note) 

P. Bonfini 1 , D. Hatzidimitriou 1 ' 3 , W. Pietsch 2 , and P. Reig 3,1 

1 University of Crete, Department of Physics, P.O. Box 2208, 71003, Heraklion, Greece 

2 Max-Planck-Institut fur extraterrestrische Physik, 85741 Garching, Germany 

3 IESL, Foundation for Research and Technology, 71110, Heraklion, Greece 



Preprint online version: July 28, 2009 



ABSTRACT 



Context. A recent survey of the Local Group spiral galaxy M 31 with XMM-Newton yielded a large number of X-ray 
sources. 

Aims. This is the second in a series of papers with the aim of identifying the optical counterparts of these X-ray sources. 
Methods. We have obtained optical spectra for 21 bright optical counterparts of 20 X-ray sources in the direction of 
M 31, using the 1.3-m Skinakas telescope in Crete, Greece. 

Results. For 17 of the 20 X-ray sources, we have identified the optical counterpart as a normal late type star (of type F 
or later) in the foreground (i.e. in the Milky Way). For two more sources there were two possible optical counterparts 
in each case, while two more objects have X-ray properties that are not compatible with the spectral characteristics of 
late type non-flaring stars. 

Key words. Galaxies: individual: M 31 - X-rays: galaxies - X-rays: stars 
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1. Introduction 

We have embarked in an effort of identifying and classifying 
the objects responsible for the X-ray emission detected in 
the XMM-Newton major-axis survey of M 31 (Pietsch et 
al. 2005, hereafter PFH2005). The PFH2005 catalogue in- 
cludes a total of 856 X-ray sources, some of which are within 
M 31 (supernova remnants - SNRs super-soft sources, X- 
ray binaries - XRBs and globular cluster sources, which 
are most likely low mass X-ray binaries within the clus- 
ter), while others are foreground stars or background active 
galactic nuclei (AGNs). A first classification of the sources 
was attempted by PFH2005 on the basis of hardness ratio 
criteria, however, several sources were not classifiable. 

The characterization of the X-ray sources involves opti- 
cal spectroscopy of candidate optical counterparts, iden- 
tified via cross-correlation with the USNO-B1 catalogue 
(Monet et al. 2003), the Local Group Survey Catalogue 
(LGS, Massey et al. 2006) and the Two Micron All Sky 
Survey (2MASS) catalogue (Skrutskie et al. 2006). The 
spectroscopic observations are conducted using different 
sized telescopes, depending on the magnitude of the candi- 
date optical counterpart. In a recent paper (Hatzidimitriou 
et al. 2006, hereafter HPM2006) we obtained optical spec- 
tra for 14 bright optical counterparts of X-ray sources in 
the direction of M 31 (along with sources in the direction 
of M33), using the 1.3-m Skinakas Telescope. The spectra 
were classified and their spectral classes compared against 
their recorded X-ray properties (fluxes and hardness ra- 
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tios) to ensure compatibility between the X-ray and optical 
characteristics. All 14 objects in the direction of M 31 were 
confirmed to be foreground stars of spectral types F0 to 
G9. 

In the present paper, we have observed another 21 
bright counterparts of 20 X-ray sources in the direction of 
M 31 from the PFH2005 catalogue, down to a limiting mag- 
nitude of Rum — 14.3 using the 1.3-m Skinakas Telescope in 
Crete (Greece). Of these X-ray sources, 18 had been classi- 
fied as foreground star candidates by PFH2005 on the basis 
of their X-ray properties and another 2 were among those 
classified as 'hard' sources (XRBs, Crab-Like SNR or AGN 
candidates). It must be noted that none of these bright 
counterparts appear in the LGS catalogue, which does not 
provide photometry for stars brighter than about V ~15 
mag. 

In a forthcoming paper (in preparation) fainter op- 
tical counterparts of X-ray sources in the direction of 
M 31, observed with the 3.5-m Telescope at Apache-Point 
Observatory, will be discussed. 

In Section 2 we describe the optical observations, in 
Section 3 the data reduction procedure, in Section 4 the 
spectral classification method and finally, in Section 5 we 
discuss the results, including individual objects of interest. 



2. Optical Observations 

The optical observations used in this study were ob- 
tained with the 1.3-m Ritchey-Cretien telescope at Skinakas 
Observatory, located on the island of Crete (Greece) and 
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were carried out during three observing runs, on September 
12-15, September 26-30, 2007 and on October 5-7, 2007. 

We used the Low Resolution Spectrograph which is in- 
corporated in the Focal Reducer instrument. A reflection 
grating of 1302 lines/mm was introduced in the collima- 
tor path, giving a nominal dispersion of 1.04A/pixel, and 
a wavelength coverage from 4680A to 676lA. The spectra 
were recorded with a 2000x800 ISA SITe CCD camera. The 
spectral region was selected to include both the H Q and Hp 
Balmer lines. 

Exposure times ranged from 1800 s to 7200 s, depending 
on the magnitude of the object and on seeing and weather 
conditions. Two or more exposures were obtained per ob- 
ject. An arc calibration exposure (NeHeAr) was recorded 
before and after each object observation. 

During the first observing run, 6 targets were observed. 
During the second run another 8 targets were observed, 
of which the candidate counterpart of source [PFH2005] 
766, USNO-B1 1314-0014843, was resolved into two stars, 
both within the X-ray source (3-sigma) error circle. During 
the final run, 6 more target spectra were obtained. Thus, 
spectra for a total of 21 candidate bright counterparts of 20 
PFH2005 sources were acquired. Table [1] provides a list of 
the optical identifications, coordinates and magnitudes of 
the objects observed along with a log of the observations. 

1. Column 1 gives the X-ray source number of the tar- 
get as it appears in the corresponding X-ray catalogue 
paper (PFH2005). As previously mentioned, for object 
766, the bright USNO-B1 star given as possible optical 
counterpart in PFH2005 is actually resolved into two 
stars on Skinakas images (with a third just outside the 
3-sigma error circle of the X-ray position). Spectra of 
both stars were obtained. Both stars are also listed in 
the 2MASS catalogue. 

2. Columns 2 & 3 give the co-ordinates of the X-ray 
sources (epoch 2000) as given in the PFH2005 catalogue. 

3. Column 4 provides the USNO-B1 identification derived 
from the cross-correlation between the X-ray catalogue 
and the USNO-B1 catalogue. 

4. Columns 5 & 6 list the B and R magnitudes of the 
corresponding USNO object (B2 and R2 in the USNO- 
Bl catalogue). The limiting magnitude of the sample of 
stars observed was R =14.3. 

5. Column 7 gives the 2MASS identification derived from 
the cross-correlation between the X-ray catalogue and 
the 2MASS catalogue. 

6. Column 8 provides the corresponding J magnitude of 
the 2MASS counterpart (as given in the 2MASS cata- 
logue). 

7. Column 9 provides the corresponding K magnitude of 
the 2MASS counterpart (as given in the 2MASS cata- 
logue). 

8. Column 10 gives the date of the Skinakas observation. 

9. Column 11 lists the total exposure time for the specific 
spectrum, while the number of exposures taken for that 
object is provided in the parenthesis following the total 
exposure time. 

10. Column 12 gives the spectral type derived in Section 4. 

A total of 14 standards of spectral types ranging from 
08 to M6 were observed during the three observing runs in 
fall 2007, with exactly the same instrument setup as the tar- 
gets. The standards were supplemented by another 15 stars 



observed with the same instrument setup in September 
2008 with spectral types between F5 and K3. These comple- 
mentary data were deemed necessary, as most target stars 
proved to be in this spectral range, and a finer spectral- 
subclass grid was needed in order to perform accurate spec- 
tral classification (see Section 4). Finally, the 14 standard 
stars from HPM2006 were also used, ranging from 09.5 to 
M6, as they were observed with the same instrument setup. 
In this manner, a grid of a total of 43 standard star spectra 
was constructed. 



3. Data reduction 

Data reduction was performed using standard procedures 
in the IRAF package 2.13-BETA2 (2006). The frames were 
bias subtracted, the sky was removed and the spectra traced 
and extracted using the all-in-one subpackage apextract. 
Arc spectra were extracted from the arc exposures, using 
exactly the same profiles as for the corresponding target 
spectra. The arc spectra were subsequently used to cali- 
brate the target spectra. Whenever more than one spectrum 
of the same object was obtained, the calibrated spectra were 
combined, after weighting them according to exposure time, 
to yield the final spectrum. 

The same reduction and calibration procedure was also 
applied to the set of the 43 spectroscopic standard stars. Q 

The signal-to- noise ratio (S/N) of the raw target spectra 
ranged between ~60 and ^200. The S/N of the standard 
star spectra ranged from ~90 to ~750. 

Figure 1 shows four examples of target spectra. 
The spectra were normalized using ordinary IRAF 
procedures. For presentation purposes, the wave- 
length region between 4800-6500 A is shown. On 
the same figure we show in red the correspond- 
ing standard star spectra of the same spectral type 
as the target spectra, shifted in y by an arbitrary 
amount. 



4. Spectral classification 

Classification of the optical spectra was achieved via cross- 
correlation against the grid of 43 standard star spectra0, 
following the procedure described in detail in HPM2006, 
i.e. each object spectrum was cross-correlated with each 
standard spectrum; the height of the corresponding cross- 
correlation peak (hereafter, CCP), which is an indicator 
of the quality of the match between the target and the 
standard star spectral features, was recorded; these CCP 
heights were plotted as a function of the spectral type of 
the standard, with the maximum of the curve yielding the 
adopted spectral type for the object spectrum. As an ex- 
ample, Figure S] shows the CCP vs. standard spectral type 
plot for object USNO-B1 1323-0017116. 



1 The spectra of the standard stars from HPM2006 were re- 
reduced following the procedure described here, for consistency, 
although no differences in the final spectra were noted. 

2 The spectral types we adopted for the stan- 
dards were taken from the Catalogue of Stellar 
Spectral Classification (Skiff, 2009) table B/mk/mktypes 
( http: / /vizier. u-strasbg.fr /viz-bin / VizieR?-source=B / mk ) ; 
when a star had more than one spectral class assigned, we 
adopted the most recent one. 
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Table 1. Log of objects observed in the direction of M 31 and derived spectral types. 
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1 X-ray source numbers from PFH2005 (M 31), i.e. the prefix [PF Z/2005] should be added in front of each number. 

2 The bright USNO-B1 star given as possible optical counterpart in PFH2005 (USNO-B1 1314-0014843) is actually resolved into 
two brighter stars on Skinakas images (and in 2MASS). 
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Fig. 2. Results of cross-correlation between the spectrum 
of object USNO-B1 1323-0017116 with each standard star 
spectrum. On the x-axis, the spectral type of the standard 
star is indicated, while on the y-axis the height of the cor- 
responding cross-correlation peak is given. 



As an independent check of the classification process, 
the spectra were also visually inspected and compared 
against the standard star spectra. 

The spectral classes derived from the CCP vs standard 
spectral type plot and from the visual inspection were in 
very good agreement in all cases. The final spectral type 
adopted for each object was the average of the two esti- 



mates (CCP and visual) and it is reported in the last col- 
umn of Table [T] 

Conservatively, the accuracy of the adopted spectral 
classes is estimated to be ~0.3 of a subclass, following the 
same arguments as in HPM2006. It is obvious that the accu- 
racy of the spectral classification will depend on the fineness 
of the grid of standard spectra used, on the signal-to-noise 
ratio of the cross-correlated spectra, as well as on spectral 
class (later spectral types are easier to classify than earlier 
types, as molecular absorption bands are very sensitive to 
temperature and thus the spectrum changes notice- 
ably between subtypes. However, these effects were taken 
into account by making the standard star grid denser wher- 
ever necessary). As in HPM2006, we tested the reliability of 
the adopted classification procedure by treating each stan- 
dard in turn as an object spectrum and by cross-correlating 
it against the rest of the standards. In Figure 2J we show 
the derived spectral type for the standards as a function of 
the reference spectral type. A least-squares linear fit yields 
a slope of 0.95 ± 0.03 (correlation coefficient R = 0.99), 
while the average disagreement between the standard and 
derived spectral type is 0.23 of a spectral type. 

In Figure 4, we plot the derived spectral types of the 
optical counterparts observed against their J — K colour as 
given in the 2MASS catalogue. The correlation is quite tight 
(correlation coefficient R = 0.94, excluding USNO- 
Bl 1314-0013202, see below), and in all but one case 
the J — K colour is in excellent agreement with the de- 
rived spectral type. The one exception is the spectrum of 
the counterpart of object [PFH2005] 273 (USNO-B1 1314- 
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Fig. 3. Derived versus standard spectral type for the stan- 
dard stars. 
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Fig. 4. Derived spectral type versus J-K colour from 
2MASS whenever available. 



0013202) which has a spectral type ^K0, while its J — K 
colour suggests type M. It must be noted that the star 
is very bright, and its recorded magnitude and colour in 
2MASS may have been affected by saturation or linearity 
effects. 



5. Discussion and Conclusions 

All spectra obtained correspond to foreground stars with 
spectral types F5 to K4. None of the stars showed emis- 
sion in H-alpha, that could be associated to flaring activ- 
ity. Figure [5] shows the frequency distribution of the de- 
rived spectral types, both from the present work and from 
HPM2006. Spectral types peak between GO and G5, with 
all but one object of type F or later. Such later type stars 
are known to be soft X-ray emitters due to coronal activity 



(e.g. Schmitt 1999, Schachter et al. 1996, Sciortino et 
al. 1993). 

In the case of [PFH2005] 766, two late type foreground 
stars of similar optical magnitude lie within the error circle 
of the X-ray source (see Table 1). Since any of these stars 
could be the source of the X-ray emission, we could not 
identify with certainty the optical counterpart. It is also 
possible that the X-ray emission results from more than 
one star, in this case. 

In the case of [PFH2005] 663, the optical counterpart 
USNO-B1 1321-0015992 (which is also recorded as a single 
star in 2MASS) is actually resolved into two close bright 
stars in the LGS catalogue. The two stars are also resolved 
on the Skinakas acquisition images, but the spectra could 
not be separated. The obtained -presumably composite- 
spectrum does not show any indication of the presence of 
two separate spectral components and was consistent with 
a K-type star. 

In order to ascertain that the optical counterparts ob- 
served are indeed responsible for the X-ray emission, com- 
patibility with the X-ray properties should be examined. 
In Table [21 we list the main X-ray properties of the sources 
(as given in PFH2005), namely the X-ray flux, f x , the log- 
arithm of the ratio of the X-ray flux to the optical flux, 
log(f x / fopt) (calculated as described in HPM2006) and the 
hardness ratios HR1, HR2, HR3 and HR4 (for the definition 
of the hardness ratios see PFH2005). 

The logarithm of the ratio of X-ray flux to opti- 
cal flux, log(f x / f pt), ranges between —3.5 and —2.8 
for the F stars in the sample, between 4.1 and —2 
for the G stars and between 3.8 and —2.4 for the K 
stars. These are consistent with the values expected 
for normal stars of types F-M (Schachter et al. 1996, 
Maccacaro et al. 1988, Vaiana et al. 1981). One of 
the G stars, the G6 star USNO-B1 1308-001239, has 
a relatively high value of log(f x / f op t) — 2 ? which 
places it at the edge of the expected range for G 
stars (Schachter et al. 1996). There is no clear corre- 
lation between log(f x / f op t) and spectral type, while 
the ranges of the possible values of log(f x / f op t) for 
the different spectral types largely overlap. Indeed, 
the soft X-ray emission caused by stellar coronal ac- 
tivity varies significantly within the same spectral 
type. This is a well known fact, which indicates that 
X ray coronal activity is the result of a combina- 
tion of several factors that determine the efficiency 
of stellar coronal heating, such as stellar magnetic 
field, stellar rotation rate, stellar age and depth of 
the convection zone (see e.g. Schachter et al. 1996). 

Figure [5] shows the hardness ratio HR2 against the X- 
ray to optical flux ratio [log(f x / fopt)], marking the differ- 
ent spectral classes in different colours. In the same figure 
we have included foreground stars in the HPM2006 sam- 
ple, excluding object [PFH2005] 464 with poorly defined 
hardness ratios (see discussion in HPM2006). All spectral 
types occupy the same locus on the log(f x / f op t)-HR2 plane. 
The objects appear to lie on a horizontal branch up to 
l°9(fx/ fopt) <-3.6; beyond this value there seems to be 
a small increase in hardness ratio. This behavior is similar 
for all the spectral types. The observed dispersion of the 
data points is expected to be affected by the unspecified 
luminosity class of the objects, as giants have on average 
higher hardness ratios and higher X-ray fluxes than dwarfs. 
Moreover high rotation can also affect the location of a 
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star in the above diagram. Two objects, [PFH2005] 40 and 
[PFH2005] 78, have hardness ratios (HR2) that appear to 
be too high for a normal late type star. Indeed, both sources 
were classified as "hard" by PFH2005. Careful inspection 
of the X-ray data did not provide any possible cause for an 
inadvertent increase of the hardness ratios recorded (such 
as close vicinity to a bright hard source), or show any X-ray 
variability. Moreover none of the stars show flaring activity 
in their optical spectra. 

In the case of [PFH2005] 40, we observed the candidate 
counterpart USNO-B1 1306-0011378. There is no other op- 
tical source within the 3-sigma error radius of the X-ray 
source position in any of the USNO-B1, 2MASS or the 
LGS catalogues. The closest neighbour is the LGS object 
J004007.76+404140.0 at a distance of 5.48" South-East of 
the X-ray position and beyond the 3 sigma error radius, at 
a magnitude of V = 21.6mag and colour B — V = 1.37. 
It is not clear that this object is a likely counterpart of 
[PFH2005] 40. There may be other sources within the 3- 
sigma error circle that are too faint to be recorded even in 
the deepest optical catalogue currently available. Figure 7 
shows the location of the X-ray source [PFH2005] 
40 along with the 3-sigma error circle, overlayed 
on a LGS i?-band image. The locations of USNO- 
Bl 1306-0011378 and J004007.76+404140.0 are also 
marked. 

In the case of [PFH2005] 78, we observed the candi- 
date counterpart USNO-B1 1304-0011697. Within the 3- 
sigma X-ray error circle, at a distance of 2.33 " South- 
West of the X-ray position, there is a V = 21.3mag 
(B — V = 0.84) object in the LGS catalogue (object 
ID: J004040.03+402549.1). Follow-up spectroscopy of this 
fainter object is necessary, in order to proceed with the 
identification of the counterpart of this sources. Figure 8 
shows the location of the X-ray source [PFH2005] 
78 along with the 3-sigma error circle, overlayed 
on a LGS i?-band image. The locations of USNO- 
Bl 1304-0011697 and J004040.03+402549.1 are also 
marked. 

In conclusion, we have obtained optical spectra for 21 
bright optical counterparts of 20 X-ray sources in the direc- 
tion of M 31, using the 1.3-m Skinakas telescope in Crete, 
Greece. For 17 of the 20 X-ray sources, we have identified 
the optical counterpart as a normal late type star (later 
than type F) in the foreground (i.e. in the Milky Way). 
For two more objects, there are two normal late type stars 
within the X-ray error circle, and therefore it was not pos- 
sible to identify with any certainty the actual counterpart 
of the X-ray emission. Two more objects have X-ray prop- 
erties that are not compatible with the spectral character- 
istics of a late type non-flaring star. In one of these cases, 
a fainter optical source within the error circle of the X-ray 
position is a possible counterpart, that needs to be con- 
firmed. 

Acknowledgements. The authors thank T. Koutentakis, N. Primak 
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M 31 


ID 1 f x (mW/m 2 ) log(f x /f opt ) 


HRl 


HR2 


HR3 


HR4 



40 


8 


.4c- 


15±le- 


•15 


-2. 


6±0. 


.1 





.12±0. 


.34 


0.56±0.18 


-0 


,24±0. 


.18 


-0 


,45±0. 


.37 


42 


1 


,0e- 


•14±le- 


•15 


-3. 


5±0. 


.1 





,53±0. 


.09 


-0.24±0.09 


-0 


.27±0. 


.13 


-0 


.59±0. 


.25 


78 


1 


,6c- 


14±2e- 


•15 


-2. 


7±0. 


.1 





.70±0. 


.18 


0.29±0.12 


-0 


,22±0. 


.12 


-0 


.35±0. 


.27 


100 


4 


,8e- 


15±lc- 


•15 


-3. 


2±0. 


.2 





.59±0. 


.18 


-0.36±0. 


.19 


-0 


,43±0. 


.30 




- 




128 


6 


,le- 


14±2e- 


•15 


-2. 


0±0. 


.1 





.56±0. 


.03 


-0.32±0. 


.03 


-0 


.84±0. 


.04 


-0 


.97±0 


.26 


129 


3 


,0e- 


15±4e- 


16 


-3. 


6±0. 


.1 





.32±0. 


.12 


-0.56±0. 


.12 


-0 


.78±0. 


.30 




- 




273 


5 


.0e- 


14±2e- 


•15 


-2. 


4±0. 


.1 





.51±0. 


.05 


-0.09±0. 


.04 


-0 


,58±0. 


.05 


-0 


.62±0. 


.21 


473 


1 


.2c- 


15±3c- 


16 


-3. 


6±0. 


.3 





.27±0. 


.23 


-0.44±0. 


.25 




- 






- 




484 


1 


.5e- 


15±4e- 


16 


-3. 


6±0. 


.3 





.04±0. 


.27 


-0.95±0. 


.20 














492 


2 


,le- 


•15±5e- 


16 


-3. 


3±0. 


.2 





.28±0. 


.14 


-0.71±0 


.27 














593 


2 


,5e- 


■15±7e- 


16 


-3. 


8±0. 


.3 





.62±0. 


.22 


-0.44±0 


.23 














611 


3 


,le- 


15±8e- 


16 


-3. 


8±0. 


.3 





.54±0. 


.21 


-0.69±0 


.23 














662 


4 


,2e- 


15±7e- 


16 


-3. 


5±0. 


.2 





.59±0. 


.13 


-0.45±0. 


.12 


-0 


.50±0. 


.24 








663 


1 


,3e- 


13±5e- 


•15 


-2. 


4±0. 


.1 





.36±0. 


.03 


-0.46±0. 


.03 


-0 


.83±0. 


.05 


-0 


.32±0 


.43 


665 


2 


,6c- 


15±8e- 


16 


-4. 


1±0. 


.3 





.51±0. 


.23 


-0.34±0. 


.23 


-0 


,83±0. 


.48 








701 


1 


,2c- 


14±le- 


•15 


-2. 


8±0. 


.1 





.42±0. 


.09 


-0.05±0. 


.08 


-0 


,63±0. 


.12 


-0 


.43±0. 


.44 


741 


5 


,2e- 


15±le- 


•15 


-3. 


3±0. 


.2 





.67±0. 


.18 


-0.31±0. 


.14 


-0 


.92±0. 


.24 








766 


3 


,2e- 


14±4e- 


•15 


-2. 


8±0. 


.1 





.53±0. 


.08 


-0.40±0. 


.08 


-0 


.59±0. 


.17 








780 


1 


,le- 


14±3e- 


•15 


-3. 


2±0. 


.3 





.72±0. 


.19 


-0.20±0. 


.22 


-0 


.60±0. 


.31 








840 


3 


.5e- 


14±le- 


•15 


-3. 


0±0. 


.1 





,70±0. 


.03 


-0.48±0. 


.03 


-0 


,90±0. 


.05 









1 X-ray source numbers from PFH2005 (M 31), i.e. the prefix [PFH2005] should be added in front of each number. 




Fig. 7. X-ray position shown with a cross (99.73% con- 
fidence error circle) of object [PFH2005] 40 overlaid on a 
R-band optical image (Massey et al. 2006). The bright stel- 
lar object with spectroscopic follow-up is USNO-B1 1306- 
0011378. The position of the Massey et al. (2006) catalogue 
object J004007.76+404140.0 is also indicated. 



Fig. 8. X-ray position (99.73% confidence error circle) 
of object [PFH2005] 78 overlaid on a V-band optical im- 
age (Massey et al. 2006). The bright stellar object with 
spectroscopic follow-up is USNO-B1 1304-0011697. The 
position of the Massey et al. (2006) catalogue object 
J004040.03+402549.1 is also indicated. 
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